Triaxial quadrupole deformation dynamics in srf-shell nuclei around ^^Mg 
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Large-amplitude dynamics of axial and triaxial quadrupole deformation in '^^'^^Mg, '^^Ne, and 
^*Si is investigated on the basis of the quadrupole collective Hamiltonian constructed with use of 
the constrained Hartree-Fock-Bogoliubov plus the local quasiparticle random-phase approximation 
method. The calculation reproduces well properties of the ground rotational bands, and P and 7 
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vibrations in Mg and Si. The 7-softness in the collective states of Mg and Ne are discussed. 
Contributions of the neutrons and protons to the transition properties are also analyzed in connection 
with the large-amplitude quadrupole dynamics. 
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PACS numbers: 21.60.-n; 21.10.Re; 21.60.Ev; 21.60.Jz 



I. INTRODUCTION 

It is known that collective deformation grows up in 
the middle of the sd-shell region. The appearance of the 
prolate ground state of ^^Mg and the oblate ground state 
of ^®Si 0-11] is associated with the shell gaps N = Z = 12 
at the prolate region and = Z = 14 at the oblate region 
in Nilsson diagram P] , respectively. Because of the shell 
gaps in the deformed regions, various shapes are expected 
to appear in the mass number region around ^^Mg and 

28Si. 

Moreover, triaxial deformation degree of freedom plays 
very important roles on the low-lying collective dynamics 
in this mass region [gj. In ^^Mg, possibility of the triaxial 
deformation in the ground states has been discussed for 
decades . The low-lying K = 2 band built on top of 
the 2 J state suggests that the triaxial degree of freedom is 
activated in the collective dynamics. In ^®Si, importance 
of triaxiality has been suggested in connection with the 
large-amplitude collective dynamics of the oblate-prolate 
shape coexistence [l^, [llj . 

In contrast to the well-developed deformations in ^^Mg 
and ^^Si, the deformation property of ^^Mg is not yet 
fully clarified. Since it is a system with iV = 14 and 
Z = 12, neutrons and protons favor different shapes 
separately. Indeed, so far many mean-field calculations 
with use of the realistic effective interactions have been 
performed for ^^Mg within an axial symmetry restric- 
tion, and they yielded a coexistence of oblate and pro- 
late shapes with an oblate minimum [l^-O- On the 
other hand, the symmetry-unrestricted mean-field calcu- 
lations using a Skyrme density functional (SkM*) [l^ 
or relativistic model show extremely triaxially soft 
potential energy surfaces. 

In the study of collective excitations in this mass 
region, the quasiparticle random-phase approximation 
(QRPA) calculations have been systematically performed 
by employing various effective interactions [Tsl . [TtI , [Tsj . 
The QRPA is a standard tool to analyze the collective 
modes of excitations. However, in order to discuss the 
low-lying collective dynamics of nuclei which are very soft 
against quadrupole deformation, one should use a micro- 



scopic theory of large- amplitude collective motion instead 
of the small-amplitude theory such as the QRPA. The 
generator coordinate method (GCM) with the restriction 
of axial symmetry [l^] , and the antisymmetrized molecu- 
lar dynamics -f multi-configuration mixing have been 
performed using the energy density functionals for mag- 
nesium isotopes and ^®Si, respectively. However, ^^Mg is 
soft against fi and 7 directions as shown in the potential 
energy surface |13 - [l5j . and therefore the triaxial degree 
of freedom in addition to the axial degree of freedom 
should be included for the description of the low-lying 
collective dynamics. 

Quite recently, the GCM calculations including axial 
and triaxial generator coordinates have been performed 
[H, [20l-|2^ for magnesium isotopes. The first applica- 
tions are concentrated on the low-lying states of ^^Mg, 
in which the small-amplitude description in the prolate 
mean field is rather good. In Ref. [la], the properties of 
the yrast states of the magnesium isotopes are discussed 
systematically. 

The quadrupole collective Hamiltonian provides a pow- 
erful theoretical tool to investigate the large-amplitude 
collective motion while taking into account the f5 and 
7 degrees of freedom [23l - [25j . Recently, on the ba- 
sis of the adiabatic self-consistent collective coordinate 
method [26l , !27| , a new microscopic method to construct 
the collective Hamiltonian has been developed, called the 
constrained Hartree-Fock-Bogoliubov plus local QRPA 
(CHFB+LQRPA) In this method, the collective 

potential is calculated by the CHFB equation, while the 
inertial functions for large- amplitude quadrupole shape 
vibration and the three-dimensional rotation are deter- 
mined from the normal modes on the CHFB state in 
(/?, 7) plane. A new point of this method is that the con- 
tributions from the time-odd mean field are taken into 
account in evaluating the vibrational and rotational iner- 
tial masses. So far this CHFB -t- LQRPA method in con- 
junction with the pairing-plus-quadrupole (P-l-Q) model 
[29I [sot including the quadrupole-pairing force has been 
successfully applied to the oblate-prolate shape coexis- 
tence in proton-rich Se and Kr isotopes [H, |3l| . 

In this paper, we analyze the role of triaxiality in con- 



2 



nection with the large-amplitude collective motion in the 
low-lying states of ^^Mg, ^^Si, ^^Mg, and ^^Ne using the 
quadrupole collective Hamiltonian calculated by use of 
the CHFB + LQRPA method with the P-^Q model. We 
also discuss the roles of neutrons and protons in N ^ Z 
nuclei on the large- amplitude collective dynamics in re- 
lation to the electric transition properties. This article is 
organized as follows. In the next section, the formulation 
of the CHFB+LQRPA method is briefly recapitulated. 
The results of the numerical calculations are presented 
in Sec. IIIIl and the role of the triaxial degree of freedom 
in these nuclei is discussed in Sec. IIVI Summary is given 
in Sec. El 

II. FORMULATION 



the Lagrange multipliers, A^^"^) = N'^'^'> —Nq'^^ are the par- 

(t) 

tide number operators measured from A^^ ^ which are the 
neutron and proton particle numbers of the nucleus. The 
operators -D2m Hermitian part of the quadrupole 

operators given by = (I>2m + £'2-m)/2. The col- 

lective potential is given by 

V{(3,^)^{^{P,^)\H\cl>{(3,^)). (6) 

On top of the CHFB state, the local normal modes are 
calculated by solving the LQRPA equations 

5(0(/3,7)|[i?CHFB(/3,7),Q^(/5,7)] 

-ip,(/3,7)|</>(/3,7))=0, (7) 



A. CHFB+LQRPA method 

The theoretical approach, the CHFB + LQRPA 
method is briefly summarized in this section. See 
Ref. (2^1 for detailed description of the method. 

The method enables us to derive the five-dimensional 
quadrupole collective Hamiltonian of the Bohr-Mottelson 
type [5, 23 25] 
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where V{(3, 7) is the collective potential in the (/?, 7) 
plane. The quantities Tvib and Trot are the vibrational 
and rotational kinetic energies. The inertial functions 
Dpp, D-yj, and Dp-y are the vibrational masses associated 
with the time-derivatives of the two quadrupole deforma- 
tion variables, /? and 7, and J7fc are the rotational mo- 
ments of inertia associated with the three components of 
the rotational angular velocities Wfc defined with respect 
to the principal axes. 

The collective potential and the inertial functions in 
the collective Hamiltonian ^ are determined microscop- 
ically in the CHFB + LQRPA method. The collective 
potential is determined by solving the CHFB equation 

(0(/3,7)|7?chfb(/3,7)I'^(/3,7)) =0, (4) 
where the CHFB Hamiltonian is given as 



(+) 

2m ' 



Tn=0,2 



(5) 



with the constraints on particle numbers and quadrupole 
deformations. Here H is the microscopic Hamiltonian, 
J, 7)) is the CHFB state, A(^)(/3,7) and Mm(/3,7) are 



-C0,-i)Q'{fin)m,i))=Q- (8) 



Here (5*(/3,7) and Pi{/3,j) are the infinitesimal genera- 
tors locally defined as functions of (/3,7). The quantity 
Ci{(3, 7) = (/?, 7) is the squared eigen frequency of the 
normal mode. We choose two collective modes from the 
LQRPA modes, following the minimal metric criterion in 
Ref. 28]- The vibrational masses Dp^^^D-yy^ and are 
determined from the transformation of the collective co- 
ordinates spanned by the two LQRPA modes into {(3, 7). 

The rotational moments of inertia are calculated by 
solving the LQRPA equations for rotation on top of the 
CHFB state. 



CHFB, *fc(/3, 7)] - -i Jk)-'ik 7)) = 0, 

(9) 



(0(/3,7)| [*fc(/3,7),4'] |0(/3,7)) =*4fc', (10) 

where ^'fc(/3,7) and Ik represent the rotational angles 
and the angular momentum operators with respect to 
the three principal axes associated with the CHFB state 
0(/3,7)}, and J^k are the LQRPA moments of inertia. 

Pauli's prescription is used to quantize the classical 
collective Hamiltonian ([T]). From the solution of the col- 
lective Schrodinger equation 



-^vib -^rc 



y(/3, 7)} *„/m(/3, 7, ^) = £^a/*a/M(/3, 7, 

(11) 



we obtain the collective wave function 5'q/m(/3, 7, i^) as 
functions of quadrupole deformations (/3, 7) and three 
Euler angles il. The collective wave function is specified 
by the angular momentum /, and its projection onto the 
z-axis of the laboratory frame, M, and a distinguishes 
the states which have the same / and M. 



3 



The collective wave function is written in the following 
form 

^^iMi^,l,n)^ ^^iK{fi,l)miMK), (12) 

-/^>0,cvcn 

where ^aiK{P, 7) is the vibrational part of the collective 
wave function, and the rotational part is written as 



{n\iMK) 



21+1 



(13) 



Here D\jj^ is the Wigner's rotation matrix and K is the 
projection of the angular momentum onto the z-axis in 
the body-fixed frame. 

The vibrational wave functions are normalized as 



d/3d7|$„,(/3,7)|2|G(/3,7)|5 = 1, (14) 



where 



|$„,(/3,7)P^ E \^o.ik{P,i)?, (15) 

-ftr >0,cvcn 

and the volume element \G{l3,^)\^dl3d'^ is given by 

|G(/3,7)|^d/3d7 = 2/3VW^(/3, l)R{P, l) sin 37^/3^7, 

(16) 

VF(/?,7) ={Dpp{p,^)D,^{p,^) - [Dp^{p,^)f}r\ 

(17) 

i?(/3,7) =i^i(/3,7)7^2(/3,7)i^3(/3,7), (18) 

where D^iP,^) are related to the moments of inertia as 
Jk (/3, 7) = 4/32i^fc (/?, 7) sin^ (7 - 27rfc/3) . 

The requantization form, the symmetries, and bound- 
ary conditions of the collective Hamiltonian are described 
in Ref. [23 1 . 

The electric properties are calculated following the dis- 
cussions in Refs. [2^ |3lj. The value of B{E2) and the 
spectroscopic quadrupole moment are given by 

B{E2-aI ^a'l') = {21 + l)-\aI\\D ^^^'^\\a' , 

(19) 

and 



QK) = \/^rj Wi^iWD'^^'^wi)- (20) 



The reduced matrix element in Eqs. and (PU)) is cal- 
culated as 

{aI\\D'^^^'^\\a'I') = I df3dj\G{f3,j)\ipifJ,,,{f3,j), 

(21) 



pLfJ'Al^,!) = V(2/ + l)(2/' + l)(-)^ 

E 



/OO^cvcn 



K k) ^"•I-K^0+ '^a'J',K' 



( f I 2 /'\ r.iE2)^ 



/ 2 



(22) 



where i a' I'if^^l) -^2 transition density. The 

quantities (/3,7) are the expectation values of the 

E2 operator in the intrinsic frame, 

i^if+'(/3,7) = ('^(/3,7)l E ei^I)i:||0(/3,7)), (23) 



where e''^ are the neutron and proton effective charges. 



n(r) 



(Dm^ + DY^j^)/2 are the neutron and proton 



and D 

parts of the quadrupole operators. 



B. Model Hamiltonian and parameters 

The pairing-plus-quadrupole model [29l [30j including 
the quadrupole-pairing interaction is adopted as a mi- 
croscopic Hamiltonian in the present work. The single- 
particle model space consists of harmonic oscillator two- 
major shells (p-shell and sd-shell) both for neutrons and 
protons. The modified oscillator values are used as the 
spherical single-particle energies ^1] . The values of the 
neutron and proton monopole pairing strengths and the 
quadrupole particle-hole interaction strengths are sum- 
marized in Table HI The strength of the quadrupole- 
pairing interaction is evaluated at the spherical CHFB 
state with use of the prescription proposed by Sakamoto 
and Kishimoto fsS^]. The interaction strengths of ^■^Mg 
are adjusted to reproduce the quadrupole deformation of 
the prolate potential minimum and the pairing energy at 
spherical shape of ^^Mg calculated with Skyrme SkM* 
density functional and the mixed surface-volume type 
pairing functional in Ref. [l8| . A simple mass number 
dependence of the interaction strengths is used to obtain 
those for ^^Mg, 24Ne, and ^^Si. Only the quadrupole 
strength for ^^Si is increased in order to adjust the de- 
formation of the oblate HFB minimum. 

The Fock term is neglected following the conven- 
tional prescription of the P+Q model. Therefore we call 
the present framework Hartree-Bogoliubov (HB) instead 
of Hartree-Fock-Bogoliubov (HFB). Following Baranger 
and Kumar [s^, the reduction factors are multiplied 
by the quadrupole matrix elements between the single- 
particle states in the upper shells, and the nuclear ra- 
dial parameter R = 1.2 A^^^ fm is used in the calcula- 
tion of the harmonic oscillator length b. In the calcula- 
tion of E2 transition strengths and quadrupole moments, 
the quadrupole operator without reduction factor is used 



[30| , and the nuclear radial parameter with the higher 
order A-dependence R = 1.12^1/3(1 + 3.84^-2/3)1/2 f^^ 
[3^ is adopted to quantitatively evaluate E2 matrix ele- 
ments. 

The two-dimensional mesh in the /3 and 7 directions is 
used to express the collective Hamiltonian in the (/?, 7) 
plane. The 60 mesh points are taken both in the range 
< /3 < /3max and 0° < 7 < 60°. As for /3„,ax, 0.6 
is used as common value, except /3max = 0.5 used for 
2^ Si, because we could not get converged solution which 
satisfies the CHB equation with four constraints at large 
deformation /3 > 0.5 near the prolate region in ^^Si. 

TABLE I: Neutron and proton monopole pairing interac- 
tion strengths gI^' and quadrupole particle-hole interaction 
strength x' = • 
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FIG. 1: (Color online) Collective potential V(/3, 7) for ^**Mg. 





Gl;'> (MeV) G^o" 


' (MeV) x' 


(MeV) 


^^Mg 


0.79 


0.83 


1.56 


26Mg 


0.73 


0.77 


1.37 




0.79 


0.83 


1.56 


28Si 


0.68 


0.71 


1.30 



III. RESULTS 
A. Collective potentials 

The collective potentials calculated for iV = Z nuclei, 
2^Mg and ^^Si, are shown in Figs. [1] and [H The collective 
potential of ^^Mg shows a prolate minimum at /? ^ 0.41, 
while that of ^^Si shows an oblate minimum at (3 ^ 0.26. 
As shown in Fig. [31 deformed shell gaps are found at 
prolately and oblately deformed regions in the present 
model. It corresponds to the appearance of the deformed 
minima in the collective potentials of ^^Mg and ^^Si. In 
2^Mg, the triaxial potential valley from the prolate mini- 
mum to the oblate region exists in the small deformed re- 
gion with 13 ^ 0.2. In ^^Si, the collective potential curve 
around the oblate minimum is steep against the triax- 
ial deformation. A prolate local minimum su ggest ed by 
other microscopic calculations for ^^Si [H, Issl |36|[ does 
not appear in the present calculation. Note that a poten- 
tial energy surface similar to the present work is reported 
in the Skyrme HF + BCS calculation [STj . 

In contrast to the deep oblate and prolate minima in 
N = Z nuclei, the collective potentials of ^^Mg and ^''Ne 
presented in Figs. H] and [S] show f3 and 7 soft situations. 
The potential minima in ^^Mg and ^"'Ne show small 
oblate deformations with (3 = 0.16 and 0.20, respectively. 
Around the potential minima, the collective potentials 
are soft against both axial and triaxial quadrupole defor- 
mations, suggesting the anharmonic situations in these 
nuclei. 

In Fig. [51 proton and neutron monopole pairing gaps 
are plotted as functions of {(3, 7). Basic characters of the 
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proton pairing gaps in ^^Mg and ^ Ne are same as that 
in ^''Mg, and those of neutron pairing gaps in ^^Mg and 
^"'Ne are same as that in ^^Si. The proton monopole 
pairing gaps in ^^Mg, ^^Mg, and ^^Ne show strong (/?, 7) 
dependence; they have the maximum at the spherical 
point, and they become zero around /? ^ 0.4 in the pro- 
late region. The neutron pairing gaps in ^^Si, ^^Mg, and 
2'*Ne becomes zero in the oblate region. These zero-gap 
regions correspond to the prolate Z = 12 and 10 shell 
gap and the oblate iV = 14 shell gap in Fig. [H respec- 
tively. Comparing the proton pairing gap for ^^Mg with 
that for ^''^Mg, and the neutron pairing gap for ^^Mg with 
that for 2®Si, an interesting feature is seen in ^^Mg; the 
zero-gap region extends widely in the (/?, 7) plane in the 
case of 2^Mg. 

In the lower panel of Fig. [31 the neutron Nilsson dia- 
gram is shown as functions of 7 at /3 = 0.3. The energies 
of the Nilsson orbits change gradually depending on 7. 
Especially, the oblate shell gap at 14 and the prolate shell 
gap at 12 more smoothly vanish in the 7 direction than 
in the /3 direction. 
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FIG. 3: Neutron Nilsson diagrams calculated within p and 
sd-shell model space as functions of /3 (a) and 7 (b). Even 
and odd parity orbits are drawn with solid and dashed lines, 
respectively. 
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FIG. 4: (Color online) Collective potential V{l3, 7) for ^"^Mg. 



B. Properties of the LQRPA modes 

As discussed in the previous subsection, one of the 
pairing features of the sd-shell nuclei is the collapse of 
the pairing gaps around the deformed shell gaps (Fig.[6|). 
In the former applications of the LQRPA equations to Se 
and Kr isotopes |28l , l3l| , where the systems are in the su- 
perconducting phase in all over the (/?, 7) region consid- 
ered, the properties and choices of the LQRPA modes in 
the normal phase has not been analyzed though they are 
interesting theoretical issues to be clarified. In Fig. [71 the 
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FIG. 5: (Color online) Collective potential V(/?, 7) for ^*Ne. 

eigen frequencies squared of the LQRPA modes Ci;^(/3, 7), 
the vibrational part of the metric W{l3,j), and the neu- 
tron and proton monopole pairing gaps Aq"'' and Aq^'' in 
^^Mg are plotted as functions of 7 along the /3 — 0.425 
line. 

Let us label the chosen two collective modes in Fig. [7] 
(a) as 'mode A' and 'mode B'. Mode A denotes the col- 
lective mode lower in energy at 7 = 0°, and this mode 
jumps around 7 = 12°, 32°, 36° and 41°. It becomes the 
second lowest mode at 7 = 60° . Mode B denotes the col- 
lective mode higher in energy at 7 = 0°, and this mode 
jumps around 7 = 8° and 25°. It becomes the lowest 
mode at 7 = 60°. 

As seen in Fig. [7] (c), in the region with 0° < 7 < 8°, 
both the neutron and proton pairing gaps vanish. In 
this region, lowest four eigen modes are close in energy 
around 1-2 MeV, and they correspond to the neutron and 
proton pairing vibrational modes at a normal phase (pair 
addition and pair annihilation) . The fifth mode is chosen 
as a collective mode (mode A). It has the 7- vibrational 
character at the axial limit, and this mode is chosen con- 
tinuously in the small 7 where the neutron is normal. 
One important character in the normal system is that 
the collectivity of the /? vibration is weakened Q . In the 
present case, the /? vibrational mode is found at around 
7 MeV. In this energy region, the /3-vibrational collec- 
tive mode is sometimes embedded in other non-collective 
modes. However, the figure shows that one can always 
find such modes using the minimal metric criterion to 
select the two collective modes. 

Around 8° < 7 < 12°, the proton becomes super- 
conducting, and the character of the low-lying modes 
changes to the neutron pair addition and annihilation, 
proton pairing vibration, and proton pairing rotation 
(Nambu-Goldstone mode). Note that in Fig. [7] (a), the 
zero-energy pairing rotational mode is not shown. 

In 12° < 7 < 25°, neutrons also become superconduct- 
ing. As the pairing gaps significantly change as functions 
of deformation in this region, the low-lying vibrational 
modes have the pairing vibrational characters, and thus 
they are not chosen to evaluate the quadrupole collective 
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FIG. 6: (Color online) Neutron and proton pairing gaps Ag 
and A^,") for ^^Mg, ^^Si, "'^Mg, and "*Ne. 

masses. In this region, the vibrational part of the metric 
increases. 

This situation changes in 25° < 7 < 41°. The pair- 
ing vibration and quadrupole vibration mix in the lowest 
LQRPA mode, and the pairing-vibrational character of 
the lowest mode decreases. Therefore, the lowest LQRPA 
mode is continuously chosen as a collective mode (mode 
B) to oblate limit. On the other hand, the other collec- 
tive mode (mode A) changes in several LQRPA modes in 
this region, because large quadrupole collectivity appears 
in these LQRPA modes. The quadrupole collectivities of 
these modes are similar, and the vibrational part of the 
metric W is continuous around the jump. 

Near the oblate region with 41° < 7 < 60°, the lowest 
two modes are chosen as the collective modes, and are 
decoupled with other LQRPA modes in energy. At the 
oblate axial limit, the lowest mode (mode B) becomes 
the 7 vibration, and the second lowest mode (mode A) 
corresponds to the (3 vibration. 

This analysis shows that the minimal metric criterion 



for selecting the two collective modes among low-lying 
LQRPA modes also works for the situation where the 
pairing gap vanishes. 



C. Collective levels 

In this subsection, we present the excitation spec- 
tra, electric transition properties, and collective wave 
functions that are obtained by solving the collective 
Schrodinger equation (fTTI) for / < 6 states. 



The excitation spectra for ^^Mg are shown in Fig. |51 
The calculation yields an yrast rotational band composed 
of 0^,2^,4:'^, and 6^ states, and an excited side band 
composed of 2^, 3]^, 4^, 5]*", and 6^ states. These two 
bands are in a very good agreement with the experi- 
mental energy levels. To analyze the structure of each 
state, the vibrational wave functions squared are dis- 
played in Fig. [9l where the /3^ factor is multiplied, which 
carries the main /3 dependence from the volume element 
|G'(/3,7)p/^d/3d7 in Eq. The members of the yrast 

band are localized around the prolate minimum, show- 
ing the prolate character of the yrast rotational band. 
The vibrational wave functions of the excited band are 
concentrated in the triaxial region around the prolate 
minimum with 7 20°. This can be interpreted as the 
7 vibrational band of the prolate yrast state. We also 
analyze the X-component fraction for these states in Ta- 
ble |TT1 The table shows that the iiT-mixing in these states 
are very small, and the results support the K = ground 
band and K — 2 excited band. These features are almost 
unchanged even with the increase of angular momentum. 

In Tables IIIII and IIVI the electric properties of the 
low-lying states are summarized. In addition to the 
theoretical results calculated with the effective charges 
(e^.g'', eg^^)=(0.5,1.5), we also list the values for pure 
neutron and proton contributions obtained by using 

{e[g\e'f^)={l,0) and (0,1). The calculated results 
reasonably reproduce the trend of the E2 transition 
strengths though they tend to underestimate the abso- 
lute values of experimental data. Using a different set 
of effective charges improves the systematical underesti- 
mation of the theoretical values. However, still there are 
some disagreement, for example, in 4^ 2^, 3^ — >■ 2^ 
and 5f 3^ transitions. One can see that the neu- 
tron and proton matrix elements are almost equal for all 
the transitions listed in the table. The prolate feature 
is shown also in the spectroscopic quadrupole moments. 
The calculated values for the yrast and the side bands 
are well described with the estimated values from the ro- 
tational collective model of a prolate deformation with 
K — and 2, respectively. This shows the rigid defor- 
mation feature in ^^Mg. 
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FIG. 7: (a) The LQRPA eigen frequencies uj^, (b) the vibra- 
tional part of the metric W, and (c) the neutron and proton 
monopole pairing gaps Aq'^' in ^''Mg plotted as a function of 7 
for a constant value /? = 0.425. In (a), the frequencies squared 
of low-lying LQRPA eigen modes are denoted with the cross 
symbols, while the two collective modes chosen at each de- 
formation to calculate the vibrational collective masses are 
denoted by the filled circles. In (b), the vibrational metrics 
are denoted by the filled circles, while the vibrational metrics 
calculated with other possible combinations of the LQRPA 
modes are denoted by the cross symbols. 



FIG. 8: Excitation spectra calculated for ^^Mg by means of 
the CHB-I-LQRPA method and experimental data 38.]. 



As mentioned in Introduction, several triaxial GCM 
calculations are performed for the study of lovif-lyina 
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states in Mg employing modern density functionals [2C 
[22| . In comparison with them, the present calculation 
gives a remarkable agreement with the experimental data 
for the ground K — band and the excited K ^ 2 band 
despite the schematic effective interaction and restricted 
model space. Especially, the agreement in the excitation 
energies are better than the GCM calculations, while that 
in the B{E2) values are worse. The small if- mixing prop- 
erties in the ground and excited bands shown in Table |TT] 
are consistent with them [2ll u3l . 



TABLE II: Probability amplitude of quantum states of ^^Mg 
for each A'-component. 



r 


K = K = 2 


K = 4 


K = 6 




1.00 










0.99 


0.01 








1.00 


0.00 


0.00 




6^ 


1.00 


0.00 


0.00 


0.00 


2t 


0.03 


0.97 






3t 




1.00 






it 


0.03 


0.95 


0.02 




5t 




1.00 


0.00 




Gt 


0.02 


0.97 


0.01 


0.00 



'St 



In addition to the ground rotational band, we show two 
rotational bands built on the 2 J and states in Fig.lTOl 
The values of B{E2) and spectroscopic quadrupole mo- 
ments are summarized in Tables |V] and IVH and the 
vibrational wave functions squared in (/3, 7) plane are 
shown in Fig. 1111 The vibrational wave functions and 
the quadrupole moments indicate the oblate deformation 
of the yrast rotational band. In Table IVIIl the K mix- 
ing in each state is listed. The yrast rotational band is 
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FIG. 9: (Color online) Vibrational wave functions squared 
/3''|$<,7(/3,7)|^ calculated for ^^Mg. 



consistent with K = oblate band. It is seen that the 
spectroscopic quadrupole moment significantly increases 
with the increase of the angular momentum. The val- 
ues of Q(4^) and Q(6^) are larger than those estimated 
from the rotational collective model using the calculated 
Q{2^) value. (In the case of the 6'^ state, 16.8 e fm^ is 
the estimated value for i^T = rotational band.) This in- 
dicates that the f3 deformation of the intrinsic state grows 
as the angular momentum increases in the yrast band as 
seen in Fig. [TTJ 

We then compare the theoretical results for the ground 
band with the experimental data. The theoretical re- 
sults overestimate both of the excitation energies and 
B{E2) values for 4^ and 6^ states. The experimental 
values of B{E2) do not follow the trend of the collec- 
tive model. This indicates that the spin alignment of 
the single-particle states plays a role in the high angular 
momentum states, and this reduces the B{E2) values. In 
the quadrupole collective Hamiltonian which we derive in 
the present calculation, the moments of inertia are eval- 
uated at zero angular momentum, and the alignment of 



TABLE III: The values of B{E2) for ^''Mg listed in units 
of e^fm^. The theoretical values are calculated with the ef- 
fective charges (e^g^e^^') = (0.5,1.5). The values calculated 
from pure neutron contribution (e^^', e^^') = (1,0) and proton 

contribution (e'fi'ie'^^) = (0,1) are also listed. Experimental 
data are taken from Ref. i3»|. 







EXP 


CHB+LQRPA neutron proton 




■>0+ 

i 


88 


63.026 


16.189 


15.614 


4^ - 


■>2+ 


160 


96.171 


24.663 


23.838 


6^ - 


■>4+ 


155 


108.032 


27.684 


26.784 


3^ - 

1 


■>2+ 

2 


239 


103.484 


26.686 


25.602 


4t - 

2 


■>3+ 




80.216 


20.674 


19.849 




■>4+ 

2 




57.085 


14.713 


14.125 


6+ - 


5+ 




47.575 


12.220 


11.786 


it - 

A 


2+ 

2 


64 


44.673 


11.448 


11.076 


K - 

i 


i 


149 


65.981 


16.944 


16.347 


6j- 


^4+ 




83.500 


21.415 


20.697 


4^- 


■>2+ 




0.011 


0.003 


0.003 


2^- 


^2+ 


15 


17.197 


4.216 


4.327 


2^- 


^0+ 


8 


4.911 


1.179 


1.244 




^4+ 




5.091 


1.241 


1.284 


3^- 


^2+ 


10 


8.180 


1.948 


2.078 




^4+ 




12.100 


2.925 


3.059 




^4+ 




0.018 


0.004 


0.005 




^2+ 


5 


3.493 


0.849 


0.882 




^6+ 




4.217 


1.018 


1.066 


5t- 


^4+ 




7.618 


1.825 


1.931 


Gt- 


^6+ 




9.756 


2.348 


2.470 




^4+ 




3.534 


0.867 


0.889 



TABLE IV: Spectroscopic quadrupole moments for Mg 
listed in units of e fm'^. See also caption in Table Hill Ex- 
perimental data are taken from Ref. [3^. 



EXP CHB+LQPRA neutron proton 



Q(2+) -16.6 


-15.7 


-7.97 


-7.81 


Qi^t) - 


-20.8 


-10.5 


-10.4 


Q(6^) - 


-23.7 


-12.0 


-11.8 


Qi^) - 


15.5 


7.87 


7.72 


Qi^t) - 











Qi4) - 


-5.85 


-3.00 


-2.90 


Qi^t) - 


-13.0 


-6.58 


-6.45 


Q{6+) - 


-UA 


-7.34 


-7.16 



the single-particle states is not taken into account. To 
include such effect by using the cranked mean field ap- 
proach is an interesting extension of the model, but is 
beyond the scope of this paper. 

The vibrational wave functions of the theoretical Oj, 
2;^, and 4;^ states have nodes in the P direction, and 
this can be interpreted as the /3-vibrational behavior on 
top of the oblate yrast states. The fact that in-band 
transitions in the ^-vibrational band B{E2;2^ — >■ 0^) 
is comparable with those of the ground rotational band 
BiE2; 2+ Of), and the positive sign of the quadrupole 
moment in 2'^ and A'^ supports this interpretation. How- 
ever, they contain K mixing and shape mixing with the 



5^ , and 6 J 



prolate region in the vibrational wave functions 

The excited band composed of 2 J, 3^, 
states are also found in the calculation. The main com- 
ponent of the vibrational wave functions is K = 2, and 
lies in the triaxially deformed region. However, there is 
no experimental information corresponding to this triax- 
ial band. Since the prolate local minimum is not found in 
the collective potential, the prolate rotational band does 
not obtained in the energy spectrum. 
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FIG. 10: Excitation spectra calculated for ^*Si by means of 
the CHB+LQRPA method and experimental data [3^ . 
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By adding two neutrons to the prolately deformed 
^^Mg, the character of collective dynamics in the low- 
lying levels drastically changes from that of ^^Mg. Figure 
1121 compares the theoretical and experimental low-lying 
energy levels in ^^Mg. In addition to the yrast band 
(0^,2^,4^ and 6^), we show a side band composed of 
2^,3^,4^,5^, and 6^ states, which are connected by 
relatively large B{E2) values (see Table IVIII| . 

The vibrational wave functions squared are shown in 
Fig. [T21 A significant difference from ^''Mg is seen in the 
deformation property of the vibrational wave functions 
of the ground 0^ state. In contrast to the well-developed 
prolate structure of the ^^Mg yrast band, the 0^ state of 
^^Mg spreads over the triaxial region from oblate to pro- 
late ones, although the shallow potential minimum is lo- 
cated at the oblate region. This indicates the very 7-soft 
character of the ground state. The members of the yrast 
band tend to localize around the prolate shape as the 
angular momentum increases. As for the excitation ener- 
gies of the yrast band, the Ex{4:f)/Ex{2f) ratio is 2.64 in 
theoretical calculation, which explains the experimental 
value 2.71 very well. As seen in Table IIXI the spectro- 
scopic quadrupole moments of the yrast band are consis- 
tent with the prolate deformation, but the absolute val- 
ues are relatively smaller than those of ^^Mg. Moreover, 
(5(6^) is almost twice of Q{2^), and this is much larger 



FIG. 11: (Color online) Vibrational wave functions squared 
/3*|$<,7(/3,7)|2 calculated for ^^Si. 



than the value of if = rotational collective model, indi- 
cating that the development of the prolate deformation 
with increase of the angular momentum. This feature is 
also seen in the ratio B{E2; 4+ 2+)/B{E2; 2+ 0+). 
The theoretical value of the ratio is 1.7, which is larger 
than the collective model value 1.43. However, the exper- 
imental value of the ratio 1.05 again indicates the effect 
of the single-particle alignment. 

The quantum states in the side band distribute widely 
in the 7 direction. This character is remarkably different 
from that of '^"'Mg, where all the members of the K ~ 2 
side band are localized in the triaxial region close to the 
prolate local minimum (Fig. [S]). In particular, an oblate 
character develops in even angular momentum states of 
the side band, and the 4j and 6j states form the two- 
peak structure in the oblate and prolate region. This 
two-peak structure indicates the 7-soft character of the 
collective potential as discussed in Ref. [soj . 
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TABLE V: The values of B{E2) for ^**Si listed in units of 
e^fm''. See also caption in Table Hill Experimental data are 
taken from Ref. fSS-l. 



FIG. 12: Excitation spectra calculated for Mg by means of 
the CHB+LQRPA method and experimental data [3^ . 
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EXP CHB+LQRPA neutron proton 




^0+ 


66.7 


41.558 


10.497 


10.354 




^2+ 


69.7 


77.816 


19.635 


19.394 




^4+ 


50.0 


98.715 


24.893 


24.608 




->2+ 


- 


65.783 


16.517 


16.422 


4- 


^3+ 


- 


41.612 


10.515 


10.366 




^4+ 


- 


27.185 


6.840 


6.782 




^5+ 


- 


30.972 


7.815 


7.719 


4- 


->2+ 


- 


59.403 


14.918 


14.828 




^3+ 


- 


62.666 


15.768 


15.633 




^4+ 


- 


98.322 


24.709 


24.538 


23+- 


^0+ 


27.8" 


44.347 


11.187 


11.053 


43+- 


■^2+ 


- 


55.761 


14.108 


13.884 


4^- 


■>2+ 


- 


0.450 


0.118 


0.111 


2+ - 


^2+ 


- 


30.834 


7.592 


7.748 


2J- 


^0+ 


- 


2.618 


0.617 


0.667 


3^- 


^4+ 


- 


8.124 


1.970 


2.051 


3^- 


->2+ 


- 


5.910 


1.387 


1.508 


4J- 


^4+ 


- 


19.928 


4.864 


5.022 


6^- 


^4+ 




0.432 


0.115 


0.106 


4J- 


-.2+ 




2.135 


0.508 


0.542 


5^- 


^6+ 




5.897 


1.408 


1.497 


5^- 


^4+ 




6.951 


1.638 


1.772 


6^- 


^6+ 




15.752 


3.813 


3.980 




->4+ 




2.836 


0.683 


0.718 


0^- 


^2+ 




56.831 


14.217 


14.205 


23+- 


^0+ 




0.078 


0.021 


0.019 


23+- 


^2+ 




2.571 


0.623 


0.649 


2^- 


^4+ 




16.159 


4.033 


4.042 



"2^ 



Go transition 




/? 0.6 ,a 0.6 



Ik AM 
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TABLE VI: Spectroscopic quadrupole moments for ^*Si listed 
in units of e fm'^. See also captions in Table Hill Experimental 
data are taken from Rcf. 38]. 

EXP CHB+LQRPA neutron proton 



Q(2+) 16 


12.0 


6.05 


5.98 


QiO - 


17.8 


8.99 


8.90 


Qiet) - 


22.2 


11.2 


11.1 


Q(2^) - 


-12.6 


-6.33 


-6.29 


0(3^) - 











Qi4) - 


-3.10 


-1.55 


-1.55 


Qi^t) - 


10.3 


5.16 


5.12 


Q(6+) - 


1.61 


0.83 


0.80 


Q(2J) - 


9.69 


4.89 


4.83 


Q(43+) - 


7.35 


3.69 


3.67 



FIG. 13: (Color online) Vibrational wave functions squared 
^^|$„7(/3,7)P calculated for ^^Mg. 



Figures[T3]and[Tn]display the energy spectra and vibra- 
tional wave functions squared for ^''Ne, and the B{E2) 
and the spectroscopic quadrupole moments are listed in 
Tables |X] and IXD The theoretical yrast band 0^, 2+, 4+, 
and (y[ has the 7-soft character similar to that of ^^Mg; 
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TABLE VII: Probability amplitude of quantum states of '^^Si 
for each _ft'2-component, where K2 is the projection of the 
angular momentum onto the symmetric axis at 7 = 60°. 





K2 = 


K2^ 


2 K2^ 


4 ATa = 6 


of 


1.00 
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2f 


0.97 


0.03 


- 


- 


4f 


0.98 


0.02 


0.00 


- 


6f 


0.99 


0.01 


0.00 


0.00 


2t 


0.11 


0.89 






St 




1.00 






4J 


0.16 


0.71 


0.14 




5t 




0.97 


0.03 




Qt 


0.15 


0.71 


0.09 


0.04 


ot 


1.00 








2f 


0.94 


0.06 






4^ 


0.67 


0.20 


0.14 





TABLE VIII: The values of B{E2) for 

~ 4 



taken from Refs. 
given in Ref. ^of 



Mg listed in units of 
lerimental values are 
3^, |40|], and the assignment of the levels are 



e^fm^. See also caption in Table Hill Experimental values are 



EXP CHB+LQRPA neutron proton 



2t- 




61.3 


52.870 


11.131 


13.953 


4f- 


-.2+ 


64.1 


90.302 


18.379 


24.070 


et- 


->4+ 




112.596 


23.016 


29.975 




■>2+ 


41.2 


74.890 


13.707 


20.568 


4f- 


■>3+ 


37.0" 


20.437 


1.366 


6.887 




->4+ 




32.073 


4.387 


9.470 


6t- 


->5+ 




19.052 


1.654 


6.156 


4j- 


■>2+ 


23.8" 


58.789 


16.303 


14.180 


5t- 


-.3+ 




70.150 


14.724 


18.530 


62+- 


■>4+ 




88.693 


23.078 


21.876 


4f- 


-+ 2j 16.0" 


2.073 


0.131 


0.704 




->2+ 


28.4 


62.940 


19.792 


14.486 


2J- 


->0+ 


1.60 


0.765 


2.040 


0.011 


3f- 


->4+ 




18.259 


5.971 


4.138 




-.2+ 


0.23" 


0.456 


3.234 


0.022 


4- 


-.4+ 




28.948 


12.298 


5.847 




■>4+ 




0.635 


0.022 


0.232 


4j- 


->2+ 




1.536 


1.755 


0.148 


6t- 


->5+ 




13.227 


4.769 


2.879 


5t- 


-.4+ 




0.773 


2.999 


0.000 


et' 


^6+ 




18.167 


9.774 


3.238 




-.4+ 




1.113 


1.006 


0.136 



TABLE IX: Spectroscopic quadrupole moment for ^^Mg listed 
in units of e fm^. See also cap tion in Table Hill Experimental 
data are taken from Ref. jSa]. 

EXP CHB+LQRPA neutron proton 



Q(2+) -13.5 


-9.07 


-2.33 


-5.27 


Q(4f) - 


-14.9 


-4.44 


-8.48 


Qi6t) - 


-17.9 


-5.61 


-10.1 


Q(2J) - 


7.17 


2.51 


3.95 


Q{3t) - 











Q{4) - 


2.74 


2.84 


0.88 


Q(5f) - 


-7.96 


-1.91 


-4.67 


Q(6+) - 


-3.25 


1.60 


-2.70 
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FIG. 14: Excitation spectra calculated for Ne by means of 
the CHB+LQRPA method and experimental data [38|. 



The two-peak structure in the 2^ , 4 J , and 6^ states 
in the side band is also seen in ^^Ne showing again the 
7-soft character as well as ^^Mg. 



The excitation energy of the 2^ state and B{E2; 2^ 



"Data taken from Ref. Uo 



of) are reproduced by the theoretical calculation. How- 
ever, the experimental energy spectrum is much more 
vibrational than the theoretical one. 



IV. DISCUSSION 



Rotational hindrance of shape mixing in Mg 
and 2*Ne 



the vibrational wave function of the Of state has a peak 
at the oblate region, but spreads over the triaxial region. 
Moreover, the gradual shape change in the yrast states 
with the increase of angular momentum is also found in 
the calculation. The oblate peak of the vibrational wave 
function in the Of state shifts to the prolate region in the 
4f state. 

The spectroscopic quadrupole moments of the yrast 
band are consistent with the prolate deformation, but 
the absolute values are smaller than -^^Mg. 



Here, we discuss the character of the vibrational wave 
functions in more detail. What is commonly seen in the 
calculated results for ^^Mg and ^''Ne is the localization 
of the vibrational wave functions squared as the increase 
of the angular momentum. While the ground Of state 
spreads over the 7 direction, the yrast band tends to lo- 
calize in the (/3, 7) plane around the prolate region, even 
though the collective potential has the shallow oblate 
minimum. 

This rotational hindrance of shape mixing is also seen 
in the cases of oblate-prolate shape coexistence [s^ |4l| 
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FIG, 15: (Color online) Vibrational wave functions squared 
/3''|$<,7(/3,7)|^ calculated for ^^Ne, 



and can be understood from the deformation dependence 
of the rotational moments of inertia. Figure [1^] shows 
the rotational moments of inertia about the intermediate 
axis, Ji(/3,7), for ^-^'^^Mg, ^'^Ne, and ^^Si, An oblate- 
prolate asymmetry is seen in the rotational moments of 
inertia; j7i(/3,7) for ^^Mg and ^^Ne becomes larger in the 
prolate side than the oblate side for the constant (3 value. 
This is the reason why the prolate shape is favored espe- 
cially for high angular momentum states. Moreover, due 
to the strong shell effect there exists a maximum point 
of the rotational moment of inertia in the {/3, 7) plane, 
which is inconsistent with the ideal irrotational moments 
of inertia proportional to the As is well known, the 
pairing correlation decreases the moment of inertia [i^ , 
and in the region where the pairing gap vanishes, the mo- 
ment of inertia becomes larger. In the case of ^^Mg, the 
neutron and proton pairing gaps vanish in the prolate 
region where the moment of inertia becomes large in this 
region (Fig, IH]). Because of the behavior of the rotational 
moments of inertia, this prolate region is favored in ro- 
tational kinetic energy, while it is unfavored in collective 



TABLE X: The values of B{E2) for ^^Ne listed in units of 
"^fm*. See also caption in Table Hill Experimental data are 



taken from Ref, [38] 



EXP CHB+LQRPA neutron proton 




^0+ 


28,0 


33,576 


14,124 


6.813 




^2+ 




56,446 


21,854 


11,906 


6^- 


^4+ 




74,510 


27,395 


16.080 


3^- 


^2+ 




45,159 


17,889 


9.426 


4.+ - 


^3+ 




6,841 


0,171 


2.579 


5^- 


^4+ 




15,843 


4,637 


3.747 




^5+ 




7,551 


0,550 


2.512 


4.+ - 


->2+ 




39,297 


15,416 


8.239 




^3+ 




44,036 


16,045 


9.540 


6^- 


^4+ 




61,481 


23,999 


12.919 



23+- 


^0+ - 


29.908 


12.458 


6.098 


4^- 


■>2+ - 


3.298 


0.419 


0.990 




-.2+ - 


53.242 


22.961 


10.675 




-.0+ - 


0.361 


0.861 


0.504 




-.4+ - 


20.227 


9.280 


3.932 


4i-- 


->3+ - 


15.732 


7.218 


3.058 


3^- 


-.2+ - 


1.487 


1.542 


1.505 


4J- 


-.4+ - 


23.575 


14.274 


3.911 


6^- 


->4+ - 


2.375 


0.325 


0.701 


it- 


-.2+ - 


0.116 


1.450 


0.030 




^► 6^ - 


16.505 


8.537 


3.008 




-.4+ - 


0.642 


2.095 


1.034 




^6+ - 


13.261 


11.282 


1.711 




-.4+ - 


0.115 


0.785 


0.005 



0.+ 



^2 

Ot 
■2| 
■4+ 



16.025 
80.124 
23.593 
0.464 
1.843 
6.076 



4.635 
23.175 
5.891 
0.044 
0.567 
2.253 



3.807 

19.034 

5.901 

0.274 

0.428 

1.306 



TABLE XI: Spectroscopic quadrupole moments for 
listed in units of e fm^. See also caption in Table Hill 



EXP CHB-f LQRPA neutron proton 



Qi^t) - 


-2.71 


1.28 


-2.23 


Q{4t) - 


-8.10 


-0.22 


-5.33 


Qiet) - 


-11.6 


-1.44 


-7.24 


Qi^t) - 


-0.21 


-1.68 


0.42 


Q(3f) - 











QiO - 


-4.68 


-1.75 


-2.54 


Qi^t) - 


-5.16 


-0.47 


-3.28 


Q(6+) - 


-9.38 


-2.80 


-5.32 


Qin) - 


-8.16 


-4.94 


-3.79 



*Ne 



potential energy, which increases as the deformation in- 
creases. As a result, the vibrational wave function for 
higher angular momentum states localizes, and the /3- or 
7-soft nature of the vibrational wave function is hindered. 

The change of the yrast state structure discussed above 
is also seen in ^*Si. In the case of ^^Si, one can see 
that the deformation of the vibrational wave function 
grows as angular momentum increases. The minimum of 
the collective potential locates around /3 = 0.26, while 
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the maximum of the moment of inertia locates around 
(3 = 0.36 in the oblate side. This is the reason why the 
deformation increases in the yrast band of ^^Si. In case 
of ^^Mg, however, the minimum of the collective poten- 
tial and the maximum of the moment of inertia coincide 
around /3 — 0.41, and such a change in the structure of 
the yrast band does not occur. 




/5 0.6 Q R 0.6 




FIG. 16: (Color online) Rotational moments of inertia about 
the intermediate axis, 



B. Analysis with mirror nucleus method 

One of the interesting issues \n N ^ Z nuclei in the 
middle of sd-shell region is the properties of the neu- 
tron and proton quadrupole transition matrix elements. 
In the mirror nucleus method [43^ . the proton matrix 
element Mp is determined from the E2 transitions as 
B{E2]I I') = M2/(2J + 1), while the neutron one 
Mn is determined from the same E2 transitions of the 
mirror nucleus as B^i„or{E2;I I') = Ml/{21 + 1). 
Although the Af„/A/p ratio should be equal to the ratio 
N/Z in the simple collective model, it has been experi- 
mentally suggested that the Mn/Mp ratio for the ground 
band transitions deviates from N/Z in some N ^ Z nu- 
clei 0, indicating possible difference between pro- 
ton and neutron shapes or shape dynamics. Moreover, 
the Mn/Mp ratio for the 2+ 0^^ transition in ^^^Mg 
is known to be extremely larger than the expected value 
N/Z in mirror nucleus method and also in the analy- 
sis of scattering reactions licl . |47| . The value of 
[Mn/Mp) /{N/Z) is 1.83 ± 0.34 [ssf. This clearly indi- 
cates the dominance of the neutron matrix element in 
this transition. In this subsection, we discuss the origin 
of the neutron dominance in the 2^ — 0]*" transition in 
26]\/jg terms of the large-amplitude triaxial shape dy- 
namics. We also discuss the yrast transition 2^ — > 0^ in 
26Mg and ^-^Ne. 



1. 2J 0+ transition m ^ Mg 

In Table IXlIl the experimental and theoretical values 
for B{E2; 2^ — > Oj^) transition are summarized. The the- 
oretical value [Mn/Mp) /{N/Z) — 2.15 reproduces the 
neutron dominance very well. Actually, the bare pro- 
ton contribution to this transition in ^^Mg is more than 
hundred times smaller than the neutron one. As seen 
from Table IIIII and |Vl such difference cannot be found in 
N — Z nuclei. The shell model value 2.10 given with 
the effective charges (e„,ep) = (0.35,1.35) also explains 
the neutron dominance of this transition. 

To analyze the mechanism of the neutron dominance 
in this transition, wc present the E2 transition density 
Pmi'W^i) defined in Eq. 1^ for ^^Mg in Fig. [HI Let 
us first compare E2 transition densities for 2^ — > 0^ 
and 2 J Of . Because of the structure of the yrast vi- 
brational wave functions, the sign of the transition den- 
sity for the former in-band transition is positive all over 
the (/3, 7) deformation, while the transition density for 
the latter transition changes its sign in the (/3, 7) plane, 
since the vibrational wave function of an excited state 
has nodes in the (/3, 7) plane. In the case of the 2^ Of 
transition as seen in Fig. [T71 the sign of the transition 
density is opposite in the prolate region and oblate re- 
gion, and this results in the cancellation after the (/3, 7) 
integration. In the case of ^^Mg, the proton matrix ele- 
ment is almost completely canceled after the integration 
of the transition density. Concerning the neutrons, the 
contribution to this transition density is relatively larger 
in the oblate region than in the prolate region, since the 
neutron favors the oblate deformation for iV = 14 system. 
This situation produces the large Mn/Mp ratio. 

This cancellation taking place in the (/3, 7) plane is the 
result of the large- amplitude collective dynamics in the 
(/3, 7) plane, especially in the 7 direction. The impor- 
tance of triaxial degree of freedom is clearly seen from 
Fig. [THl where the E2 transition density is plotted as a 
function of 7 for a constant value of /3. It is seen that the 
proton transition density is almost anti-symmetric with 
respect to 7 = 30°, while the asymmetry is present for 
neutron transition density. 



TABLE XII: The values of B{E2;2+ 0+) for ^^Mg and 
^*Si listed in units of e^im^. Theoretical values for '^^Mg are 
calculated with the effective charges (e^g' , e^g' ) = (0.5,1.5), 
while those for '^^Si are calculated assuming the mirror sym- 
metry. The bare neutron and proton contributions calculated 
with (e^g\e^g) — (1,0) and (0,1) are also listed. Experimen- 
tal data and shell model calculation are taken from Refs. [ssl ], 
and ':48j, respectively. 

EXP CHB-HLQRPA neutron proton SM" 
^"Mg 1.60±0.32 0.765 2.040 0.011 3.28 

^^Si 7.32±2.28 4.82 0.011 2.040 19.7 
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FIG. 17: (Color online) E2 transition density ^*p^^^',^, (/3, 7) for 2+ 0+ and 2+ 0+ transitions in ^'^Mg. Three sets of 



effective charges (e^J^', e^^') =(0.5,1.5), (1,0) and (0,1) are used for left, middle and right panels, respectively. 
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FIG. 18: (Color online) E2 transition density P'^p^L'i'iPn) 
for 2^ of transition in ^'^Mg at ^ = 0.345 is plotted as 
function of 7. The red solid, green dashed, and blue dotted 
lines show the calculations with three sets of effective charges, 
(eeff ,e^ff (0.5,1.5), (1,0), and (0,1), respectively. 



proton deformation is expected in the ground state of 
24Ne m. 

We evaluate the M„/Mp ratio by using the mirror nu- 
cleus method. In TableKftll B{E2; 2+ 0+) values for 
^^Mg, ^^Si, ^^Ne, and ^^Si are summarized. The relative 
magnitudes of E2 transition probability for mirror pairs 
cannot be satisfactory reproduced by the theoretical cal- 
culation both for A — 24 and 26 systems. 

The theoretical value of the ratio {Mn/Mp)/{N/Z) for 
26]\/[g jg 0.81. The calculated value is smaller than unity, 
and qualitatively reproduces the tendency of the experi- 
mental value. The shell model value [l^ gives smaller ra- 
tio 0.69 than the present calculation. For ^^Ne, the ratio 
is calculated to be 0.86. The results for ^^Ne fail to quan- 
titatively describe the suppression of the ratio extracted 
from the central values of the experimental B{E2). For 
more detailed discussions, precise measurements of the 
E2 transition strengths for ^''Ne and ^^Si are required. 



2. 2J Q{ transition in ^ Mg and Ne 

In contrast to the large {Mn / Mp) / {N / Z) for the 2^ 
0^ transition discussed above, the experimental value of 
the ratio {Mn/Mp)/{N/Z) for the 2'l 0^ transition 
in ^^Mg is 0.92 ± 0.05, which is close to unity. This in- 
dicates that a simple collective model picture with the 
usual assumption that the radius and the deformation 
for neutrons are consistent with those for protons is ex- 
pected to hold. In the case of ^^Ne, the experimental 
value of this ratio is 0.59 ± 0.11. The possible suppres- 
sion of the ratio from unity may suggest that the simple 
picture does not hold for neutrons and protons in this 
system, and the smaller neutron deformation than the 



TABLE XIII: The values of B{E2;2+ -s- 0+) for ^''Mg, ^^Si, 
^^Ne, and '^^Si are summarized in units of e^fm*. Theoretical 
values for ^^Si and ^''Si are calculated assuming the mirror 
symmetry. See also caption in Table IXlIl Experimental values 
and shell model values are taken from Refs. 38, 44, 49], and 
Ref. [3], respectively. 



EXP 


CHB+LQRPA neutron proton 


SM 


^"Mg 61.4± 1.8 


52.870 


11.131 13.953 


57.1 


^''Si 70.5± 6.9 


47.226 


13.953 11.131 


36.5 


^*Ne 28.0± 6.6 


33.576 


14.124 6.813 




2^Si 19.1± 5.9 


48.197 


6.813 14.124 
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V. SUMMARY 

Large-amplitude triaxial quadrupole deformation dy- 
namics in the low-lying states of srf-shcU nuclei, ^''Mg, 
^*Si, ^^Mg, and ^^Ne are analyzed on the basis of the 
quadrupole collective Hamiltonian derived microscopi- 
cally from the CHFB + LQRPA method. 

As for the N ~ Z systems, the calculation reproduces 
the prolate rotational band and the 7 vibrational band 
in ^''Mg, and the oblate rotational band and the /3 vi- 
brational band in ^^Si. As fov N ^ Z systems, ^^Mg 
and ^^Ne, the collective potentials are shown to be soft 
against the (3 and 7 deformations, and the large shape- 
fluctuation in the 7) plane is found in the vibrational 
wave functions of the ground states. 

The yrast bands show rotational hindrance of the 
shape mixing, and the states localize around the pro- 
late region as the angular momentum increases. The 
neutron and proton quadrupole matrix elements are ana- 
lyzed for N ^ Z systems. The neutron dominance in the 
2^ 0^ transition in ^^Mg is explained in terms of the 
large-amplitude collective dynamics in the 7-direction. 
The neutron and proton matrix elements for 2f — > 0^ 
yrast transition are analyzed with use of the mirror nu- 



cleus method for ^^Mg and ^^Nc, Also in other N Z 
nuclei, differences in the behavior of neutrons and pro- 
tons in large-amplitude shape dynamics are expected to 
be interesting. 
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